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ABSTRACT 

We study lower and upper bounds on the parameters for stochastic 
state vector reduction, focusing on the mass-proportional continuous spon- 
taneous localization (CSL) model. We show that the assumption that the 
state vector is reduced when a latent image is formed, in photography or 
etched track detection, requires a CSL reduction rate parameter A that is 
larger than conventionally assumed by a factor of roughly 2 x 10^^^, for a 
correlation length rc of 10~^cm. We reanalyze existing upper bounds on the 
reduction rate and conclude that all are compatible with such an increase in 
A. The best bounds that we have obtained come from a consideration of heat- 
ing of the intergalactic medium (IGM) , which shows that A can be at most 
~ 10^=*=^ times as large as the standard CSL value, again for rc — 10~^cm. 
(For both the lower and upper bounds, quoted errors are not purely statis- 
tical errors, but rather are estimates reflecting modeling uncertainties.) We 
discuss modifications in our analysis corresponding to a larger value of rc- 
With a substantially enlarged rate parameter, CSL effects may be within 
range of experimental detection (or refutation) with current technologies. 
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1. Introduction 



Stochastic modifications of the Schrodinger equation have been intensively studied as 
models for objective state vector reduction [1]. As currently formulated, the reduction rate 
parameters for these models are many decades smaller than current experimental bounds, 
and will not be detectable in planned nanomechanical and gravitational wave detector exper- 
iments [2]. In order to motivate further experimental searches for stochastic modifications of 
Schrodinger dynamics, it is important to have lower bounds on the stochastic model param- 
eters, below which one can assert that known measurement processes will not occur. Setting 
new, more stringent lower bounds, while at the same time reanalyzing and improving upper 
bounds, is the aim of this paper. 

Conventional lower bounds on the stochastic parameters are based on an idealized 
measurement model in which the experimenter reads out results from the position of a macro- 
scopic pointer. In such measurements, the detection and amplification processes needed to 
get a pointer readout are necessarily linked. Our focus in this paper is on a different type 
of experiment, in which first a latent image is formed, either in a photographic emulsion or 
a solid state track detector. Only long after latent image formation is amplification brought 
into play, in the form of development of the photographic plate, or etching of the track de- 
tector. A qualitative discussion of latent photographic image formation was given in 1993 by 
Gisin and Percival [3], who consider detection to have occurred already at the microscopic 
level, before the amplification associated with development. To quote them, "Of particular 
importance is the study of the formation of the latent image in photography, for this is 
not only the most common quantum detection technology, but it also shows unequivocally 
that amplification up to the macroscopic level is quite unnecessary for the formation of a 
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permanent classical record of a quantum event, by contrast with the example of the pointer, 
which is so often used." In this paper we shall put the discussion of Gisin and Percival on a 
quantitative footing, within the framework of the continuous spontaneous localization (CSL) 
model with mass-proportional couplings. 

This paper is organized as follows. In Sec. 2 we give a quick review of the CSL 
model with mass-proportional couplings, and in particular give rate formulas needed for the 
subsequent discussion. In Sec. 3 we discuss the formation of latent photographic images 
in the Mott-Gurney model, and show that with standard parameter values, the CSL model 
predicts a reduction rate that is a factor of order ~ 2 x 10^^^ too slow, as compared with 
the estimated rate of latent image formation. In Sec. 4, we make analogous (but cruder) 
estimates for solid state etched track detectors, and again conclude that the usual CSL 
parameter values cannot account for latent image formation. 

In Sec. 5, we examine whether various upper bounds on the CSL parameters allow 
a substantial enlargement in the reduction rate parameter. Processes considered include 
FuUerene diffraction, supercurrent persistence, proton decay, spontaneous radiation from 
germanium, cosmic IGM heating effects, and planetary heating. Our conclusion in all cases 
not involving heating (in some cases disagreeing with previous analyses) is that an increase 
of A by a factor of 10^^ is allowed by experimental data. Heating of the IGM places a 
stronger constraint on A, allowing an increase by a factor of ~ 10^"^^ over the standard value. 
Although planetary heating nominally places a much more stringent bound, we argue that 
competition with molecular collision effects, which are 28 orders of magnitude larger and are 
strongly dissipative, invalidates this bound. 

In Sec. 6, we discuss modifications in our analysis resulting from changing the cor- 
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relation function g{x), and from changing the value of the correlation length rc from the 
standard value 10~^ cm assumed in CSL analyses. In Sec. 7 we discuss implications of our 
suggested new values for the CSL parameters for experiments to directly test for CSL ef- 
fects, focusing on large molecule diffraction, superconductor current decay, nanomechanical 
and gravitational wave detection experiments, and the CoUett-Pearle [4] proposal to observe 
rotational Brownian diffusion. We also briefly consider the competition of molecular col- 
lisions with CSL effects. Finally, in Sec. 8 we discuss and summarize our results, and in 
particular, we note that our key assumption, that latent image formation already constitutes 
measurement, is subject to direct experimental test. 

2. The mass-proportional CSL model 

We begin by stating some standard formulas of the mass-proportional CSL model, 
drawing heavily on the review of Bassi and Ghirardi [1]. The basic stochastic differential 
equation of the model is 



d\m) 



-^Hdt + J d^x{M{x) - {M{x)))dB{x) j d^x{M{x) - {M{x))fdt 



\m) , 
(1) 

with dB{x) a Brownian motion obeying 

dtdB{x) = , dB{x)dB{y) = -fd^{x - y)dt , (2) 

with (M) denoting the expectation of M in the state \il^{t)), and with the operator M(x) 
given by 

M{x) = m],' J d^yg{x - y) Y^^sNsiy) . (3) 



In Eq. (3) the sum extends over particle species s of mass rris and with number density 
operator Ns{y), while uin is the mass of the nucleon and g{x) is a spatial correlation function 
conventionally chosen as 

9i^) = e-^^/'>\ jd'xg{x)^l . (4) 

(In Sec. 6 we will show that changing the functional form of the correlation function does 
not substantially alter our conclusions.) This correlation function can be interpreted as the 
functional "square root" of the correlation function for a noise variable that couples locally 
to the mass density X^s That is, writing dC{y) — J (fxg{x — y)dB{x), which has 

the spatial correlation function dC{y)dC{z) = jdt J d^xg{x — y)g{x — z), the noise term in 
Eq. (1) takes the locally coupled form J d^ydC{y)ni^Ylis^s{Ns{y) — {^s{y)))- Thus, the 
basic parameters of the model are the strength of the Brownian process 7, and the correlation 
function width parameter a. 

It is convenient to introduce two further parameters that are defined in terms of 7 
and a. Since a has the dimensions of inverse squared length, we define a correlation length 
rc (sometimes denoted by a in the CSL literature) by 

rc = («)-'/' , (5a) 

which is conventionally assumed to take the value of rc = 10^^ cm. Additionally, we 
introduce a rate parameter A defined by 

A = 7(^)'^' = 7/(87r^/VS) , (56) 

which is conventionally assumed to take the value A = 2.2 x 10^^^ s^^, giving 7 the value 
10-30cm3s-i. 
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We can now state two key rate formulas that we will need for the subsequent analysis. 
According to Eq. (8.15) of Bassi and Ghirardi [1], the off-diagonal coordinate space density 
matrix element (^|p|0) for a single nucleon approaches zero exponentially with a reduction 
rate Tr given by 

r« = A(l - e-^V(44)) , (6a) 

which for i comparable to or larger than rc, can be approximated as 

Tr^X . {6b) 

For n nucleons within a radius smaller than the correlation length, this rate is multiplied 
by n^; for N groups of nucleons separated by more than the correlation length, this rate is 
multiplied by N, and for particles of mass nip the rate is multiplied by {mp/mN)^, giving 

TR^Xri'Nimp/mNy ■ (6c) 

The formulas of Eqs. (6a-c) will be the basis of our reduction time estimates for latent image 
formation. Before proceeding further, we note that there are other natural definitions of a 
reduction rate; for example, if it is defined by the rate of approach to zero of the variance in 
the coordinate, as in Adler [1], rather than the rate of vanishing of the off-diagonal density 
matrix element, then Fr is double that given by Eqs. (6a-c). 

A second important formula gives the rate of secular center-of-mass energy gain, as 
a result of the Brownian process, for a body comprised of a group of particles of total mass 
M. This is given by the formula (Bassi and Ghirardi [1], Pearle and Squires [5], Adler [6]) 

dE _?> M 

This formula will be used to set upper bounds on the reduction rate parameter A. 
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3. Latent image formation in photography 

The case of latent image formation in photography is of particular interest because 
the photographic process has been intensively studied, both theoretically and experimentally 
Survey accounts of what is known are given in the books of Mott and Gurney [7] and Avan 
et al [8], and in the review articles of Berg [9] and of Hamilton and Urbach [10]. 

A photographic emulsion consists of grains of silver halide suspended in gelatine. 
The typical grain size is between 1/10 of a micron and a few microns, or in other words, 
from 10~^ cm to a few times 10"'' cm in diameter. Grains are typically spaced around 1 to 
2.5 microns apart in the gelatine. The basic microscopic theory originated by Gurney and 
Mott [11] is illustrated graphically on p 75 of ref [8], for a typical AgBr emulsion. The steps 
as envisaged by Gurney and Mott are as follows. First a photon is absorbed in the grain, 
giving rise to an electron and a hole. The electron gets trapped on the surface of the grain, 
and the hole produces a neutral Br which is absorbed on the surface. An interstitial ion of 
Ag"^ then diffuses to the trapped electron and forms a neutral silver atom. A second photon 
is absorbed, giving rise to another electron and hole; the hole produces a second neutral Br 
which joins with the first to produce a Br2 molecule, which eventually diffuses out of the 
grain into the gelatine, while the electron combines with the neutral Ag on the surface to 
give an ion Ag~. (As shown in Fig. 5.4 of Hamilton and Urbach [10], the bromine that 
has diffused into the gelatine typically moves up to around a micron away from the grain; 
this will figure in our discussion of Sec. 6.) Finally, a second interstitial Ag"*" diffuses to 
the surface to join the Ag~, forming a molecule Ag2. This process is then repeated until 
a cluster of free silver atoms (a silver speck) has been formed on the surface of the grain, 
making it developable. As discussed in detail in ref [8], an alternative microscopic model has 
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been proposed by Mitchell, in which crystalline imperfections play an important role, and in 
which the order of steps differs from that in the theory of Gurney and Mott. However, such 
details are not relevant to the estimates that we shall make, which only use the number of 
ions that move and the distance that they travel, but not the order in which the motions 
take place, details of trapping, etc. 

There seems to be general agreement that around 3 to 6 silver atoms are needed to 
give a grain a 50% probability of being developable, and about 30 silver atoms are needed 
for a grain to be certain to be developable. In using grains to define a particle track in a 
100 micron width emulsion, a few grains developable by random processes are encountered, 
so we shall assume that a minimum number of about N — 20 developable grains is required 
to define a track. Assuming that 30 Ag and Br atoms (with respective atomic weights of 
108 and 80) move a distance greater than or of order rc in the ionic step of latent image 
formation, we have for the number of nucleons that move n = 30(108 + 80) = 5640. So 
from Eq. (6c), we have for the reduction rate Fr for the process of latent image formation 
producing the track, 

Tr = Xn^N = 2.2 X lO'^V^ x 5640^ x 20 = 1.3 x lO'V^ . (8) 

On the other hand, Mott and Gurney [7] estimate in their Table 28 the rate of 
growth of a silver speck. Scaling their room temperature example to a single speck of radius 
0.2 X 10~^cm, which fits comfortably within a grain of diameter 10~^cm, one finds a rate for 
accumulating 30 silver atoms of around 30s~^. This rate should be considered uncertain by 
two orders of magnitude in either direction; for example, using the ion step time of 1/25, 000s 
discussed on p 248 of Mott and Gurney, the rate to accumulate 30 silver atoms would be 
about lOOOs"^, while applying the results of Table 28 to a speck of radius 0.2 x 10~^cm, 
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corresponding to just a few silver atoms, would give a rate to accumulate 30 atoms of 0.3s~^. 
Note that since the accumulation of silver atoms on a speck is a sequential process, the time 
for building up the speck can be much longer than the duration of the exposure that initially 
produces electrons and holes in the interior of the grain. 

If we identify the time to form a developable speck with the total time to form the 
latent image, we see that a reduction rate of 3 x 10^='=^s~^ is required for reduction to be 
completed during latent image formation. Thus, based on this estimate, the rate of Eq. (8) 
is too small by a factor of 2 x 10^=*=^; in other words, the CSL rate parameter A would have to 
be increased by this factor to account for latent image formation (for rc fixed at 10~^cm). 

Let us now address a possible objection to this conclusion, that we have considered 
only the mass motion associated with registration of the latent image, but have neglected 
possible motions of environmental particles induced by the registration process. The first 
thing to be said is that photographic emulsions work perfectly well in vacuo, so it suffices to 
consider just the motions of particles within the emulsion. When a photon is absorbed by an 
atom, the atom will recoil, potentially altering the phonon distribution, and thus displacing 
all the atoms in the detector. For a photon of energy a few eV absorbed by a silver hahde 
molecule, this recoil will have an energy of order 10~^°eV, which is much less than the 
typical Debye temperature energy Ed ~ lO^^eV. Thus the parameter Er/Ed governing the 
probability of a phonon- free absorption is of order 10~^, and so the recoil will be absorbed, 
with probability 1 — O(10~^), by a translation of the whole detector, without the emission 
of phonons. Hence for relevant elapsed times t, the recoil distance I will be much smaller 
than the correlation length, and we can use the expansion of Eq. (6a) given in Eq. (9) below. 
Given a total detector mass of NnniN, the recoil distance will he £ — pt/{Nn'mN), with p 
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the incident photon momentum, and so the reduction rate associated with detector recoil 
will be Tr = [X/{4:N)]\pt/{mNrc)]'^- Notice that the factor has canceled out, and the 
rate is proportional to the inverse of the emulsion size, as reflected in the factor N~^. For 
a 3eV incident photon, a correlation length rc = 10~^cm, and an elapsed time t of l/30s, 
this estimate gives ~ 0.5A^~^ x 10~^s~^ for the standard CSL value of A. For a minimal 
detector with = 20 we obtain a reduction rate about a factor of two larger than that 
estimated in Eq. (8), but for a one cubic centimeter emulsion, with ~ 10^^ nucleons, the 
effective N in the recoil estimate is 10^^, and the recoil induced reduction rate is neghgible. 

Once the incident photon has been absorbed, the remaining steps in latent image 
formation involve transport and diffusion of electrons, holes, or silver and bromine ions. 
The electron liberated from a silver halide molecule by the absorbed photon will initially 
have a velocity larger than the mean thermal velocity, and will emit phonons until it slows 
down to thermal velocity This effect is calculated in Appendix D, where we estimate that 
it leads to a reduction rate F^ ~ 10~^^s~^, nine orders of magnitude smaller than the result 
of Eq. (8). Assuming that interstitial ions are already thermalized, and that they have a 
dispersion relation of the form E — /2m* with m* an effective mass, the corresponding 
phonon emission reduction rates due to ionic motions are much smaller, because only the 
tail of the ionic thermal momentum distribution is kinematically allowed to emit phonons. 
Once thermal equilibrium is established, only the difference between the thermal motions in 
superposed states in the wave function can contribute to reduction. This difference should 
be smaller than twice the thermal motion in one of the states of the superposition, and an 
estimate given in Appendix D shows that this bound on the thermal contribution to the 
reduction rate is similar in magnitude to the result of Eq. (8). Finally, the motion of a 
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silver (or bromine) ion from the interior of a grain to the surface requires a compensating 
back motion of the whole emulsion if no momentum is transferred to the surroundings. For 
example, if a silver ion diffuses a distance rc, to keep the center of mass of the emulsion 
fixed, it must translate by 108rc/{Nn), leading to a reduction rate Tr ~ x 10~^^s~^ 
for the standard CSL value of A. This is neghgibly small even for = 1, that is, for an 
emulsion consisting of a single grain. In sum, these estimates suggest that collective motions 
of the atoms in the emulsion, induced by the process of latent image formation, do not alter 
our conclusion that the standard CSL value of A is not large enough to account for latent 
image formation. 

4. Latent image formation in etched track detectors 

In this section we turn to the consideration of a second type of latent image formation, 
in which etchable defects are formed in solid state detectors. Using related parameters, 
we will also estimate the reduction time associated with electronic motions resulting from 
passage of the charged particle through the detector. 

Detailed accounts of etched track detectors are given in the books of Durrani and 
Bull [12] and Avan et al [8], from which we draw the needed parameters. Passage of a charged 
particle results in a cylindrical region of lattice distortion around the charged particle path, 
with a typical radius of 6 to 12 nm (i.e., 60 to 120 A, or 0.6 to 1.2 xlO~^ cm), which we 
take as giving the displacement i in Eq. (6a). Since this i is considerably smaller than rc, 
we expand the exponential to get the associated reduction rate 

rR^Xn'Nf/{irl) , (9) 

where as before, n is the number of nucleons within a correlation length, and N is the number 
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of correlation-length groups of nucleons in the particle track. Let us estimate the length of 
a short ionic track length to be 10 microns, giving N = 100. For n, we take the number 
of nucleons in a cylinder of diameter i — 10"^ cm, and length 10~^ cm, assuming a density 
of 10^"^ nucleons per cubic centimeter. This estimate gives n ~ 3 x 10^. Thus, with the 
standard value A = 2.2 x 10~^^s~^, we get a reduction rate of ~ 7 x 10~^s~^. Various 
numbers for the time t for a track to form are given in the two books cited. A minimum 
for t is given by the lattice vibration time of 10~^^s; according to Table VIILS of Avan et 
al, establishment of thermal equilibrium requires a time in the range of 10~^^s to 10~^s, 
and establishing chemical equilibrium requires a time of a few times 10^® s to much longer, 
depending on the material. Thus, these time estimates indicate that the standard value of 
A gives a reduction rate that is too small by a factor of 10^^ to 10^"^, if thermal equilibrium 
is taken as the criterion for latent image formation, and a factor of around 5 x 10^ or less, 
if chemical equilibrium is used as the criterion. 

Let us next address a potential objection to this conclusion, and to the similar one 
reached in Sec. 3 in the case of photographic emulsion grains, that we have neglected the 
mass transport associated with the ionization produced by passage of the charged particle. 
However, the initial ionization consists almost entirely of the motion of electrons, and ac- 
cording to Katz and Kobetich [13], the bulk of the energy deposited within 20 A of the track 
comes from 5-rays with energy less than O.lkeV. As one might expect, these electrons have 
ranges of order 40A or less, so the geometry of electron displacements is essentially the same 
as the geometry of nucleonic displacements. Since the electron-associated reduction rate is 
smaller than that calculated above by a factor (mg/mAr)^, the electrons that are initially 
ionized in the course of track formation give a negligible correction to the estimate made 
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above. 



5. Upper bounds on the CSL pcirameter 

In the previous two sections, we have concluded that if state vector reduction occurs 
when a latent image is formed, the reduction rate parameter A must be much larger than 
conventionally assumed. We turn now to an examination of whether such an enlarged A is 
allowed by various experimental constraints. We continue to assume for this discussion that 
the correlation length rc is 10~^cm, an assumption to which we shall return in Sec. 6. 

5A. Fullerene diffraction experiments 

We begin our survey of experimental constraints with the fullerene diffraction exper- 
iments of Arndt et al [14], Nairz, Arndt and Zeilinger [15], and Nairz et al [16], which we 
previously discussed in Sec. 6.5 of Adler [1]. In these experiments, molecules with n ~ 10^ 
nucleons are observed to produce interference fringes when diffracted through gratings with 
spacings of 1 to 2.5 x 10~^cm, of order the correlation length rc- Prom Eq. (6c), we get a 
reduction rate with the standard A of = 2.2 x IQ-^^s"^ x 1000^ ~ 2 x IQ-^^s'^ Since 
the beam transit time in these experiments is of order 10~^s, the interference fringes will not 
be spoiled provided that Tji < 10^s~^, permitting the rate parameter A to be up to 5 x 10^^ 
times as large as the standard value. 

5B. Decay of supercurrents 

The decay of supercurrents in the mass-proportional CSL model has been calculated 
by Buffa, Nicrosini and Rimini [17], as also reported in Rimini [18], extending to the CSL 
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case earlier work by Rae [19]. Equation (5.9) of ref [17] gives a fractional supercurrent decay 
rate, using the standard CSL parameters, of 5.1 x 10~^^s~^. Since the current experimental 
limit is [20] about 10~^year~^ = 3 x 10~^^s~^, the rate parameter A is permitted to be up to 
10^^ times as large as the standard value. (This limit will be increased when recombination 
processes leading to the formation of Cooper pairs are taken into account, since these are 
not included in the calculations of refs. [17], [18], and [19].) 

The review of Bassi and Ghirardi [1] quotes in detail only the numbers obtained by 
Rae, although referring to the later work of refs [17] and [18]. Rae assumes a fractional 
decay rate given by A, for which he takes the value 10~^^s~^, as opposed to the CSL value 
of 2.2 X 10~^^s~^. The CSL model calculation in refs [17] and [18] decreases the estimate 
given by A by a factor of {rckp)'^ — 0.6 x 10~^ (with kp the Fermi momentum) arising 
from indistinguishability of electrons, and a factor (me/Mjv)^ — 0.3 x 10~®, arising from the 
reduced stochastic coupling of the electron in the mass-proportional scheme. Assuming that 
each spontaneous localization breaks a Cooper pair, and that Cooper pairs are not recreated, 
this gives an overall supercurrent decay rate of 2.2 x lO^^^s"^ x 0.2 x 10~^ = 4.4 x 10~^'''s~^ 
essentially the same number given in ref [17]. 

The indistinguishability factor {rckp)'^ can be understood heuristically from the 
fact that localization within a radius rc involves a momentum transfer 5k ~ l/r^, and so is 
forbidden by the exclusion principle except for electrons within 5k oi the Fermi surface, which 
are a fraction of order 5k/ kp of all electrons. Since for temperatures well below the critical 
temperature all electrons in the Fermi sea contribute to the supercurrent, and not just those 
near the Fermi surface, spontaneous localization can thus affect only a fraction {rckp)^^ of 
the electrons in the supercurrent, and hence this factor appears in the supercurrent decay 
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rate. 

5C. Excitation of Bound Atomic and Nuclear Systems 

Pearle and Squires [5] have estimated tlie rate of bound state excitation in the mass- 
proportional CSL model and applied it to various physical systems. Their Eq. (12) gives for 
the rate P of internal excitation of atoms 

P ^\{me/mNf{ao/TcY , (10) 

with oo the atomic radius ~ 10~^cm, which with the standard value of A gives P ~ 0.7 x 
10~^^s~^. Requiring that cosmic hydrogen not dissociate during the hfetime of the universe 
of 4 X lO^^s gives an upper bound on A of 4 x 10^'' times the standard value. 

Also in ref [5] , Pearle and Squires apply an analogous formula to proton decay, with 
ao now a nuclear radius of 10~^^ cm, and with me replaced by the quark mass rrig. Using 
the constituent quark mass rUg ~ mAr/3, this gives P ~ 10~^°s~^, while using current quark 
masses of around 10 MeV gives P ~ 10~^^s~^. Requiring P to be smaller than the proton 
decay rate, which is known to be less than 10~^^year^^ = 0.3 x 10~^°s~^, allows A to be 
up to ~ 10^ times the standard value if constituent quark masses are used in the estimate, 
and up to ~ 10^^ times the standard value if current quark masses are used. However, 
these limits are overly restrictive because the formula of Pearle and Squires assumes that no 
selection rules are at work. This is true in the case of atomic dissociation, but not in the 
case of decay of a proton. Because free quarks are confined, and because fermion number 
is conserved modulo 2, allowed proton decay modes must contain at least one lepton, which 
is not originally present as a constituent of the proton. Hence the proton decay rate, which 
involves a wave function overlap squared, must be further suppressed by (m^v/A)^, with A 
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the minimum scale at which physics beyond the standard model is found. Making the very 
conservative estimate A ~ 250GeV (the electroweak scale) gives an additional suppression 
of 3 X 10~^°; hence proton decay allows A to be up to at least 10^^ times the standard value, 
using the constituent quark mass estimate. 

Yet another application of the excitation rate formula is given by CoUett and Pearle 
[4], to an experiment in which the rate of spontaneous 11 keV photon emission from ger- 
manium, as monitored in 1 keV bins, has been bounded by 0.05 pulses/(keV kg day). The 
assumption here is that nuclear excitation will have a high probability of knocking out a 
Is electron, with resulting emission of an llkeV photon as the remaining electrons cascade 
downward to fill the vacant orbit. To get an estimate we use Eq. (10) with mg replaced by 
mjv, and Oq taken as the nuclear radius given by 1.4 x 1Q~^^A^^^ cm, with A the atomic 
number. Using ^4 ~ 73 for germanium, which has 8.3 x 10^^ atoms/kg, and equating P to 
the rate limit for the first IkeV bin, gives the bound A < 6 x 10~^s~^, which is ~ 3 x 10^^ 
times as large as the standard CSL value of A. (This bound is a factor of 10^ higher than the 
one quoted by CoUett and Pearle [4]; where they get A~^r^ > 2 x 10~^^cm^s, our evaluation 
from their numbers gives A~^r^ > 2 x lO^^^cm'^s.) Yet another analysis of nuclear excita- 
tion has recently been given by Pearle [4] , who used the CSL model to estimate the rate of 
quadrupole excitation of a germanium nucleus to its first excited state. Pearle shows that 
for this process, the experimental data on photon emission require A to be less than ~ 10^^ 
times its standard value. 

5D. Radiation by Free Electrons 

Pu [21] has calculated the rate of radiation of photons of energy k by free electrons 



17 



that is induced by the stochastic term in the Schrodinger equation. Including the factor 
{me/niNy called for in the mass-proportional coupling case, his result becomes 



An alternative derivation of this result is given in Appendix A, obtained by calculating 
the mean squared acceleration produced by the stochastic term, and then using this in the 
standard formula for the power radiated by an accelerated charge. In evaluating his result 
numerically, Fu takes = 1/137.04, whereas the standard Feynman rules that he uses 
require the identification eV(47r) = 1/137.04. He also uses the value A = 10 ^^s ^, so to 
correct for his evaluation of and to correspond to the CSL value of A, we multiply the rate 
result of his Eq. (4.1) by 47r/4.5 = 2.8. Comparing with the experimental bound on photon 
emission from germanium quoted above, this gives A < 1.7 x 10~^^s~^, which is only ~ 10^ 
times as large as the standard CSL value of A. 

However, in the special case of mass-proportional coupling that we are considering 
here, this bound is significantly modified. Although the valence electrons in germanium are 
quasi-free, they are still charge neutralized by the core consisting of the inner electrons and 
nucleus. The core has a center of mass motion with a characteristic lattice vibration period 
of order 2 x 10~^^s, much longer than the 4 x 10~^^s period of an emitted 11 keV photon, 
so the core motion should give at most an order 10~^ correction to the rate of 11 kilovolt 
photon emission by the core induced by the stochastic term in the Schrodinger equation. 
In the mass-proportional coupling CSL model, since the core radius is much less than the 
correlation length, the core can be treated as a point particle with the stochastic term acting 
only on its center of mass. Moreover, the acceleration of the core induced by the stochastic 
term will be the same as that of electrons (the factor relating acceleration to force is 



dr{k) 

dk 



e^A 



(11) 
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canceled by the m factor in the couphng to the noise), and so the core motion will tend 
to cancel the radiation field produced by the electron motions. The cancellation in the far- 
zone radiation field is incomplete because of two types of corrections. The first are order 
O'o/rc corrections, refiecting the variation of the correlation function over the germanium 
atom radius Oq ~ 10~^cm, which gives a correction of order (oo/rc)^ ~ 10~^ times Fu's 
estimate in the radiated power. The second are retardation corrections of order v/c, with 
~ 0.3 X 10~^c the thermal velocity of the valence electrons, which gives a correction of 
order v'^/c^ ~ 10~^ times Fu's estimate in the radiated power. (There will also be a cross- 
term correction in the radiated power of order {ao/rc){v/c) ~ 10~^ times Fu's value.) Thus, 
taking the effect of charge neutralization into account, Fu's calculation shows that A can be 
at most ~ 10^^ times as large as the standard CSL value of A. 

5E. Heating of Protons 

We consider next the heating of protons by the stochastic noise, as given by Eq. (7). 
Taking M = rriN, 'Vfe get from Eq. (7) and the standard CSL parameters, 

^ = -A^^ = 6.8xlO-26eVs-i = l.lxlO-3Vgs-^ . (12) 

Hence over a ten billion year period (3.15 x 10^ ''s), ignoring dissipation, a proton would 
gain an energy from this effect of ~ 2 x 10~^ eV, corresponding to a temperature of ~ 
2 X 10~^ degrees Kelvin. If A were larger that the standard CSL value by a factor of 2 x 10^, 
the minimum value of the range inferred from our latent image formation analysis, the 
temperature increase would be ~ 7 x 10^ degrees Kelvin over the age of the universe. 

We use this number to get a number of different bounds. The first is obtained by 
supposing that if A is much larger than the standard value, all of the energy gained by 
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protons in the universe is thermalized into low energy photons. The total energy per unit 
volume released this way must not exceed a small fraction, say 0.1, of the energy per unit 
volume in the 3 degree (~ 2.6 x 10""^ eV) microwave background radiation. Taking account 
of the fact that there are ~ 10^ microwave background radiation photons per proton in the 
universe, this gives the bound that A must be less than ~ 10^^ of the standard CSL value. 

A more precise way to get a cosmological bound from proton heating is to utilize the 
fact that there have been detailed experimental and theoretical studies of low density regions 
of intergalactic space containing the intergalactic medium (IGM) [22]. The IGM consists of 
highly ionized hydrogen, with a typical temperature of around 2x10^ degrees Kelvin [22] as 
measured between redshifts oiz — A and z — 2. The IGM is heated by radiation from various 
astrophysical sources (stars, supernovas, quasi-stellar objects) and is cooled by adiabatic 
expansion of the universe, and by recombination cooling of the plasma. We can get an upper 
bound on the proton heating rate, by ignoring possible astrophysical heating mechanisms, 
and assuming that all of the heating in fact comes from the CSL heating effect, which must 
be equated to the cooling rate rate from adiabatic expansion and recombination cooling to 
explain the observed temperature equilibrium. For the highly ionized IGM, recombination 
cooling is less important (by a factor of around 6 at ;s = 3) than adiabatic expansion, so to 
to get an estimate we shall ignore the influence of recombination cooling, and just consider 
the effect of adiabatic cooling as follows. Under adiabatic expansion, the temperature varies 
with redshift z as T oc (1 + z)^, and so the thermal energy loss of a proton obeys dE/dt — 
{d/dt){3/2)kT = {3kT/{l+z)) \dz/dt\. At z=3 (for Hubble constant Ho = 71 km/s/Mpc and 
for matter and dark energy fractions Qm = 0.26 , Qa = 0.74) one has \dt/ dz\ = 0.8 x lO^years. 
Assuming a temperature T at 2; = 3 of 2 x 10^ degrees Kelvin, this gives an energy loss rate 
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of 0.5 X lO^^^eVs"^. Equating this to an energy gain from proton stochastic heating, from 
Eq. (12), gives a value of A equal to 8 x 10*^ times the standard CSL value as an upper bound. 

A second IGM bound can be obtained from measurements at z ~ 0, which are less 
precise than those in the range 2 < z < i. Ricotti, Gnedin, and ShuU [22] give a z = 0.06 
mean temperature of ~ iQ^ '^^o.b (jggrggg Kelvin (with one sigma errors), which limits A to be 
at most i0''-2='=°-^ times as large as the standard CSL value. Combining the z — 3 IGM bound 
obtained in the preceding paragraph with this one, and rounding the error in the exponent 
to an integer, we get an overall IGM bound on A of ~ 10^^^ times the standard CSL value. 

Yet another cosmological bound comes from considering interstellar dust grains. Here 
we compare the rate of energy accumulation, given by Eq. (12), with the rate of energy 
radiation by the grain. Assuming 10^"^ nucleons per cubic centimeter, the rate of volume 
energy production for the standard CSL value of A is ~ 7 x 10~^eVs~^cm~^. Radiation from 
dust grains does not follow the Stefan-Boltzmann law (power oc area x T^), but rather scales 
with the dust grain volume and the fifth power of the temperature [23], according to the 
following formula giving W, the rate of energy radiated per unit volume of grain, 

iy = 327r24.94?^/^' . (13) 

Here h and c are the Planck constant and the velocity of light, Tg is the grain temperature, 
and k' is the imaginary part of the refractive index. Taking typical values Tg ~ 20 degrees 
Kelvin and ac' ~ 0.05, this formula gives Vl^ ~ 2 x 10^^eVs~^cm~^. Hence the dust grain 
energy balance implies only a weak bound, that A can be at most of order 10^^ times its 
standard CSL value. 

We finally consider bounds on energy production obtained from the energy balance 
associated with planetary heat fiows. Table 6.3 of de Pater and Lissauer [24] gives the lumi- 

21 



nosity to mass ratio L/M for solar system objects. The lowest ratios in the table, < 4, 4, 6.4 
in units of 10~^ergg~^s~^, come respectively from Uranus, carbonaceous chondrites, and 
Earth. Comparing with the nucleon heating rate of Eq. (12), which is ~ 7 x lO'^^ergg^^s"^, 
and ignoring the fact that the nucleons in earth are not isolated, we would conclude that A 
can be at most 5 x 10^ times the standard CSL value. However, this bound is dubious, for the 
following reason. In a recent paper, Bassi, Ippoliti and Vacchini [25] , following up on earlier 
work of Halliwell and Zoupas [26] and Gallis and Fleming [27] (see also Benatti et al [28]), 
point out that in stochastic models with a dissipative form, the rate of heating of nucleons 
is not constant, but approaches zero after a finite time, so that the energy gained through 
the heating effect approaches a finite limit as opposed to increasing indefinitely. Even when 
the fundamental stochastic reduction process is non-dissipative, for nucleons in the earth 
the reduction process is a much weaker effect (by about 28 orders of magnitude, for A at 
the IGM upper limit) than ordinary coUisional decoherence [29], which is dissipative. Hence 
any stochastic heating effect should have equilibrated to zero very rapidly as a result of the 
influence of molecular collision effects, and so Eq. (12) cannot be compared to the planetary 
heat flow to set a bound on A. 



6. Possible modifications in our analysis 



We discuss in this section possible modifications in our analysis. We consider first 
a modification that might eliminate the upper bound obtained from IGM heating. The 
heating bound could be eliminated if, as discussed in ref [25], the stochastic Schrodinger 
equation itself were modified to include dissipative effects. If the energy gain by a proton 
from stochastic effects were to saturate, independent of the value of A, at a level significantly 
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below 10^ degrees Kelvin, then IGM heating would not lead to a bound on A. However, as 
emphasized in ref [25], for a dissipative stochastic model to be a serious competitor to the 
standard CSL model, a way has to be found to give dissipative reduction a field theoretic 
formulation, so that Fermi and Bose statistics for identical particles can be properly taken 
into account, as is done in the CSL model. This is an important open problem needing 
further study, but in the absence of its solution we continue to assume that the stochastic 
Schrodinger equation is exactly nondissipative in form. 

We consider next the effect on the heating bound of changing the correlation function 
g{x) and the correlation length rc- The first comment to be made is that changing the form 
of g{x), as long as it remains a function only of the magnitude of x, has no qualitative effect 
on the analysis. For example, suppose that instead of a Gaussian, one takes the simple 
exponential form 

^(x) = e-l^l/'-^/(87rrS) , (14a) 

where the normalizing factors have been chosen to ensure that J d^xg{x) = 1. At first sight, 
one might think that such a choice would lead to the appearance of single inverse powers of 
rc in various expressions, such as the small ^ expansion of Tr and the stochastic heating rate 
(see Eqs. (9) and (7) respectively). This, however, is not the case, because as is made clear 
in Eq. (8.14) of the review of Bassi and Ghirardi [1], the function g[x) enters the analysis 
only through the integral 

G{y) = j d'xg{y - x)g{-x) . (146) 

When g(x) is inversion invariant {g{x) = g{—x)), the order y term in the Taylor expansion 
of g{y — x) vanishes on integration over x, and thus G{y)/G{0) begins with an order y^ term. 
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For example, the choice g{x) of Eq. (14a) leads (with s — \y\/rc) to 

C.(.)=^(l + . + .V3)e- 

(14c) 

= — [l-sV6 + sV24-sV45 + 0(s^)] , 
D47r 

in which an odd power of \y\ appears first at fifth order! Hence with alternative choices of 
the correlation function, the qualitative form of Eqs. (9) and (7) remains unchanged, apart 
from a rescaling of rc by a numerical constant. See Weber [1] for a further discussion of 
alternative choices of the correlation function. 

Continuing, then, with the Gaussian choice for g{x), let us consider the effect of 
changing rc from the conventional, and rather arbitrary choice, of rc — 10~^cm. Clearly, if 
rc were increased to rc ~ 10~^cm = l/i, the IGM heating rate given by Eq. (7) would be 
reduced by two orders of magnitude, and the upper bound on A would be raised to 10^°^^ 
times the standard CSL value. 

We next analyze the effect on our latent image bounds of increasing rc- We consider 
first the case of formation of a photographic latent image. For a small AgBr grain of diameter 
i = 10~^cm, motion of the Ag and Br ions within the grain would be suppressed, according to 
Eq. (9), by a factor |(-^/rc)^ = 1/400. However, since the bromine atoms that leave the grain 
typically move a distance of up to a micron into the gelatine, the suppression factor for these 
atoms would be only of order | dxx^ — 1/12. So with n = 30 x 80 = 2400, and as before 
= 20 and A = 2.2 x 10-^^s-^ our revised estimate from Eq. (9) is = 2.2 x 10-^°s-\ 
requiring A to be 1.4 x 10^^^^ times larger than the standard CSL value. 

Turning next to latent image formation in etched track detectors, when rc is in- 
creased by a factor of 10, the value of n used in Eq. (9) also increases by a factor of 10, but 
N decreases by the same factor. Therefore the combination n^N/r'c is a factor of 10 smaller 
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than estimated previously, and so an increase of A by a factor of < 5 x 10^" is required for 
the reduction rate to equal the latent image formation rate, using attainment of chemical 
equilibrium as the criterion. 

To conclude this section, we note that the physiology of the eye gives a possible 
argument suggesting that rc should not be larger than 1 or 2 microns. Consider a dense 
array of photoreceptors, which one is trying to optimize both with respect to rapidity of 
response, which we assume to be given by the reduction rate Tr, and the ability to resolve 
fine detail. If the diameter of each photoreceptor is £, there will clearly be a tradeoff between 
the two desired attributes of the array, as we vary £. If i is much smaller than rc, the 
response rate of each detector element is suppressed by a factor of with a factor of i"^ 
coming from the exphcit i'^ in Eq. (9), and an additional factor of coming because the factor 
n in Eq. (9) scales as £. (Since we are considering here the rate for a single photodetector 
element, the factor N in Eq. (9) will be 1). On the other hand, the spatial resolution of 
the array varies as Hence one gains more in response time than one loses in spatial 
resolution as £ is increased, until £ reaches rc- Once £ is larger than rc, the increase in the 
response time becomes only linear in £ rather than quartic (c.f. Eqs. (6a) and (6c), with n 
now a constant and N increasing in proportion to £.) Thus increasing i beyond rc produces 
a gain in response time that is less than the corresponding loss of spatial resolving power. 
So these considerations suggest that an optimized photoreceptor array should have £ ~ rc- 

Assuming that the evolution of the human eye has optimized it with respect to both 
response time and resolving power, and using the empirical fact [30] that the diameter of 
rods in the retina is 2/i, this reasoning suggests that the correlation length rc is unlikely 
to be larger than of the order of a micron. However, some caveats are in order. First, 
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since physiological constraints prevent eucaryotic cells from being much smaller than of the 
order of microns, this argument cannot be used to disfavor smaller values of rc, such as the 
conventionally assumed rc — 0.1 micron. Also, in applying the optimized array argument 
to the rods of the eye, we have implicitly assumed that with the enhanced reduction rate 
suggested by latent image formation, reduction occurs directly within the individual rods, 
rather than requiring molecular motions in the optic nerve as well. As discussed briefly 
in Appendix C, estimates based on the known amplification chain in the rods support this 
assumption, but are incompatible with reduction occurring in the conformational change of 
an individual rhodopsin molecule. However, these estimates (if not overly optimistic) suggest 
such a high reduction rate in the rod amplification chain that reduction rate is no longer a 
relevant factor in determining the optimized array dimensions, which then undermines the 
argument just given relating the diameter of rods to the value of rc- 

7. Implications of a larger reduction rate 

We turn now to a discussion of experimental implications of a greatly enhanced 
value of A. We divide our analysis into two parts, first considering those experiments that 
do not rely on the secular increase of energy associated with the CSL process, and then 
considering experiments that depend on this secular energy increase (in which case the 
analysis possibly would be modified in dissipative extensions of the CSL model.) We consider 
two possible choices of enhanced CSL parameters, (1) rc — 10~^cm, A = 4 x 10~^°s~^, and 
(II) rc — 10~^cm, A = 3 X 10~^s~^. These correspond respectively to the latent image lower 
bounds of 2 X 10^ and 1.4 x 10^ times the standard CSL value, and are near the lower end 
of the respective ranges permitted by the IGM heating upper bounds. 
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7A. Experiments not based on seculeir energy increase 

Fullerene diffraction For case (I), with a grating of lO^^cm — re, setting = \v? — 
10^s~^ gives a value of n ~ 5 x 10^ at which washing out of the diffraction pattern would set in. 
For case (II), with a grating oi I — 2.5 x 10~^cm = .25rc, setting F^^ = \n^i^ j (4r^) = lO^s"-^ 
also gives a value of n ~ 5 x 10^ at which the diffraction pattern would start to wash out. 
Thus, diffraction experiments with projectiles of molecular weight of order 500,000 (a factor 
of 500 beyond what has been achieved so far) would confront the CSL model, when the 
parameter values are minimally enhanced to account for reduction in latent image formation. 

Supercurrent decay For case (I), from the calculation of refs [17] and [18] we find a 
supercurrent decay rate of 10~^^s~^ ~ 1/(3 x lO^^years), while for case (II), we find a 
supercurrent decay rate of 8 x 10~^^s~^ ~ 1/(4 x 10^°years). A direct measurement of the 
supercurrent decay time constant (as opposed to a time constant inferred from an improved 
resistivity measurement based on a short measurement time) would be needed to confront 
the CSL model with enhanced parameter values. Again, we note that these estimates ignore 
recombination processes in which Cooper pairs are created. 

Mirror deflection experiment The papers of ref [2] give a detailed analysis of a mirror 
deflection experiment proposed by Marshall et al [31], as expressed in terms of the param- 
eter 77 governing the small displacement CSL equation. For the geometry of the proposed 
experiment, with a cubical mirror of side S and density when S » tq the parameter 77 
is given by the formula 

77 = %'Kr%\S'^D^ . (15) 

When S is of order rc, this formula has corrections of order unity, but we shall continue to 
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use it to give a rough order of magnitude estimate. For case (I) we find tfiat 77 is a factor of 
2 X 10^ larger than in ref [2], giving a fringe visibihty damping factor e~^, with A ~ 0.04. 
The thermal decoherence background has been estimated recently by Bernad, Diosi and 
Geszti [32] as At ~ 0.5, so the case (I) CSL effect is not detectable. For case (II), we find 
that ?7 is a factor of 1.5 x 10^^ bigger than than in ref [2], corresponding to A ~ 3 x 10^, 
indicating complete suppression of the interference fringes, a large and detectable effect. 
Thus, even a version of the Marshall et al experiment that is a factor 100 times less sensitive 
than envisaged in their proposal would be of considerable interest, provided that thermal 
decoherence backgrounds can be kept small. 



We discuss here experiments that utilize the fact that any non-dissipative Brownian 
process, including the CSL stochastic noise, will produce a secular increase in system energy 
and a corresponding increase in the rms values of translational and rotational displacements. 
Should it be possible to formulate a satisfactory dissipative version of CSL, the predictions 
for these experiments would be altered when the observational time exceeds a characteristic 
time for dissipative effects to set in (see ref [25]). Having stated this caveat, we proceed to 
give the results expected from the standard version of CSL, with enhanced parameters as 
suggested by latent image formation. 

The CoUett— Peairle rotational diffusion proposal CoUett and Pearle [4] have proposed 
searching for the mean square rotational diffusion A^csl of a suspended disk of order rc in 
dimensions. We write the result of their Eq. (C.6) in the form 



7B. Experiments utilizing secular energy increase 




(16a) 
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with /rot a dimensionless function of the ratios of disk dimensions to rc given in ref [4], 
with / the disk moment of inertia, and with the "standard quantum hmit" A6'sql given by 



For readers not famihar with standard quantum hmits, Eq. (16b), and also the analogous 
result for the standard quantum limit on measurement of the translational coordinate of 



Eq. (16b) is substituted into Eq. (16a), which is why it does not appear in Eq. (C.6) of ref 
[4].) Following CoUett and Pearle, we assume a disk radius L = Ire and a disk thickness 
h — 0.5rc, for which /rot — 1/3, and the moment of inertia is given by / = Ml? j^^ 
with M the disk mass. Assuming a disk density of lO^^cm"^, we then find, for case (I), 
A^csl/A^sql = 6.6 X lO^t, and for case (II) A^csl/A^sql = 1-8 x lO^t, with t in seconds 
in both formulas. For estimates of the rotational diffusion produced by Brownian motion, 
which places severe pressure and temperature constraints on the environment of the disk, 
see Eq. (6.4) of ref [4]. Since the Brownian rotational diffusion scales as the product of 
the square root of the pressure and the fourth root of the temperature, increasing the CSL 
effect by a factor of order 10^ to 10^ makes the vacuum and cryogenic constraints on the 
experiment much less severe. 

Nanomechanical oscillator We consider next the nanomechanical resonator reported by 
LaHaye et al [33], and discussed from the viewpoint of CSL effects in Adler [6], to which 
the reader is referred for details. Recalculating the numbers found in ref [6] to reflect an 
enhanced CSL parameter A (and in case (II), an increased rc), we find that the increase in 
the harmonic oscillator occupation number N over a time equal to the inverse of the noise 
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(166) 



a free mass, are given in Appendix B. (Note that the moment of inertia / cancels when 



bandwidth is of order 10~^ in case (I), and of order 10~^ in case (II). Similarly, the root 
mean square deviation of the quantum nondemolition variables Xi^2 over a time equal to the 
inverse of the noise bandwidth is of order 3 x 10~^ of the corresponding standard quantum 
limit in case (I), and of order the 0.3 times the standard quantum limit in case (II). Hence 
the enhanced CSL effect would still be not accessible in case (I), and would come close to 
being detectable in case (II). 

Gravitational wave detectors LIGO and LISA We consider finally gravitational wave 
detector experiments, also discussed from the viewpoint of CSL effects in ref [6] . The Ad- 
vanced LIGO Interferometers [34] are expected to approach the standard quantum limit in 
position measurement accuracy, while the Laser Interferometer Space Antenna (LISA) [35] 
should achieve a position accuracy of around 10^ times the standard quantum limit. Recal- 
culating the numbers from ref [6], we find that for LIGO the root mean square stochastic 
deviation in the test mass coordinate over a time interval of 1/70 s will be ~ 0.02 times the 
standard quantum limit in case (I), and ~ 1.4 times the standard quantum limit in case 
(II). For LISA, the corresponding figures over a time interval of 10^ s are ~ 6 x 10^ times 
the standard quantum limit (that is ~ 6 times the expected position accuracy) in case (I), 
and ~ 5 X 10^ times the standard quantum hmit (~ 5 x 10^ times the expected position 
accuracy) in case (II). Thus, the enhanced CSL parameter values give an effect that may be 
barely visible in LIGO in case (II), and should be readily visible in LISA in both cases (1) 
and (II). 

The background from molecular collisions In designing experiments to look for CSL 

effects connected with the secular increase in energy and associated coordinate deviations, 
one will have to make sure that similar effects associated with molecular collisions are smaller 
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in magnitude. There are two ways of estimating such effects. The first is to use the decoher- 
ence calculation of Joos and Zeh [29] (with later corrections [29] in the numerical coefficient), 
which gives a rate parameter A for decrease in the off-diagonal density matrix element in a 
coordinate basis, that is analogous to the CSL parameter ri/2. The second is to use standard 
classical Brownian motion formulas, as is done in the paper of CoUett and Pearle [4]. In 
Adler [36] we calculate and compare the mean square coordinate fluctuation predicted by 
the two methods, for the case of isotropic hard sphere scattering. We show that the two 
calculations give identical results, and so either method can be used to estimate collision 
backgrounds that may mask CSL effects. 

8. Summary and discussion 

To summarize, we have found that if latent image formation constitutes measure- 
ment, then the parameters of the CSL model for objective state vector reduction must be 
much larger than conventionally assumed. The enhanced parameter values needed are com- 
patible with empirical upper bounds, and suggest that CSL effects may be within reach of 
experimental detection within the next decade or two. 

Further experiments on latent image formation would also be of great interest, to 
firm up (or falsify) assumptions that we have made in our analysis. It would be very use- 
ful to have more accurate information on the time needed for a latent image to form, in 
both photography and etched track detectors, since estimates of this time (or rate) were a 
crucial input for our analysis. Our central assumption, that latent image formation, and 
not subsequent development or etching, constitutes measurement, should also be subject to 
experimental test. In principle, a photographic emulsion or an etched track detector could 
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be used as a "which path" detector in one arm of a quantum interferometer; if latent image 
formation constitutes measurement, then the interference fringes should be destroyed by the 
presence of the emulsion or detector through formation of the latent image, which is a "fos- 
silized record" of the path through the interferometer. According to the view expressed by 
Gisin and Percival [3] and adopted in this article, collapse onto a definite outcome should 
occur once the latent image has formed, and should not depend on a subsequent stage of 
amplification by development or etching. 
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Appendix A: Radiation by a free charge 

We give here an alternative calculation of the radiation by a free charge, obtained by 
using the basic dipole radiation formula P = ^e^a^, with e^/ (47r) = 1/137.04 and with a the 
acceleration of the particle of charge e. For a time average over the stochastic process, is to 
be interpreted as E[{x)\ with E[ ] the stochastic expectation. Rather than proceeding from 
the CSL equations, we use a modification of the CSL equations that has a unitary dynamics 
but leads to the same density matrix evolution; this substitution is allowed because the 
expectation £■[ ] of quantum transition probabilities is expressible directly in terms of the 
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density matrix. (This substitution is used by Fu [21] as the basis of his calculation.) The 
advantage of doing this is that the modified dynamics has a Heisenberg picture formulation 
with a self-adjoint Hamiltonian, which permits a direct calculation of x from the Heisenberg 
equations of motion. 

Specifically, corresponding to the CSL model we have the effective coordinate repre- 
sentation Hamiltonian 

H-^-JtJ dh{m/mN)dB{z)g{z - x) , (Al) 



from which we calculate 



X ^ - Ux,H]^ — , 
n m 



i h f . 

X = - -\p/m,H] = ^^^^ / d^zdB{z)da:g{z - x) , 

with dx a vector gradient. Thus for E[{x)'^] we get 



= 4^ / d'zMz - x)r = -^J^ . (A3) 



Using 1/dt = 5{0) = {2tt) ^ dk = tt ^ dk, and substituting into the dipole radiation 
formula, we get 

47r2 mj^rl; Jo 

in agreement with Eq. (11) when specialized to units with h — 1. The total radiated power is 
infinite, but since CSL is a nonrelativistic theory, the calculation is reliable only for /c << m. 
With a cutoff taken at A; = m, the radiated power is smaller than the power gain through 
stochastic heating (see Eq. (7)) by a factor e^/(37r^). 

We also note that exactly the same result is obtained from the so called QMUPL 
model, which corresponds to the leading small-displacement Taylor expansion of the CSL 
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model (and of the original Ghirardi-Rimini-Weber [37] model as well), as discussed and 
further referenced in refs [2]. In this case the effective Hamiltonian is 

h 

with {dBf = dt and r] = A/(2r^), and one immediately finds 

tun dt 

Since the QMUPL model of Eq. (A5) is one-dimensional, one must multiply by a factor of 3 
in forming the expectation of the squared acceleration from Eq. ( A6) , giving the same result 
as obtained in Eq. (A3) from the CSL calculation. As one would expect, the integral over 
the squared gradient oi g{z — x) in Eq. (A3) is just the one encountered in evaluating r] from 
the second order expansion of the function G{y) of Eq. (14b), in the course of deriving the 
QMUPL model from the small displacement expansion of the CSL model. 

Aside from giving a simpler derivation of Fu's result (which he obtained by using 
Feynman rules to calculate the S matrix for the radiation process) , the derivation given here 
shows that the radiation is directly attributable to the noise-induced acceleration. Moreover, 
the radiation has a coefficient that is independent of the particle mass m, because the 1/m in 
the relation between x and p cancels against the m coming from the noise coefficient m/miq 
included in the mass-proportional coupling scheme. This cancellation of m dependence is why 
the radiation field from the germanium core largely cancels that from the valence electrons. 

Appendix B: Standeird quantum limits 

We derive here the standard quantum limit for measurement of the angle 9 through 
which a suspended disk rotates about the axis of suspension, taken here as the z axis. We 

34 



then state by analogy the corresponding standard quantum hmit for the measurement of a 
translational coordinate. 

Let Lz — —ihd/d9 be the component of angular momentum around the z axis, which 
obeys the commutation relation 

[e,L,] = zh . {Bi) 

We shall only consider very small angular displacements, and so we neglect complications 
associated with the 2tt periodicity of 9 and the corresponding integer quantization of L^. 
Prom the commutation relation of Eq. (Bl), one gets an uncertainty relation of the usual 
form, 

AML,>| . {B2) 

Suppose at t — a. measurement of 9 is made to accuracy A^. This induces an uncertainty 
ALz > h/{2A9), and since AL^ = lAuj^, with / the moment of inertia of the disk around 
the z axis and uj^ the corresponding angular velocity, the squared uncertainty in the angular 
coordinate of the disk at time t is 

{A9{t)f ^{A9f + {tAuj,f ^{A9f + {thf/{2IA9f . (S3) 

Minimizing this expression with respect to A^, we get the standard quantum limit on the 
measurement accuracy of 9 at time 

{A9{t)Y > {A9sql)\ A9sQL^{ht/iy/' . (S4) 

In an entirely analogous fashion, from the commutation relation [x, p] = ih, together 
with p — Mv with M the mass and v the velocity, one derives the standard quantum limit 
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for measurement of a translational coordinate, Axsql = {Ut/My/'^. For further details of 
standard quantum limits, see ref [38]. 

Appendix C: Where does reduction occur in the visual system? 

We give here some estimates for the implications of an enhanced A parameter for 
the human visual system. Attention was first drawn to the issue of possible state vector 
reduction in the visual system by Albert and Vaidman [39], with a response in Aicardi et 
al [39], as reviewed in Bassi and Ghirardi [1]. An issue raised in this early work is: where 
does one expect reduction to occur in the visual system - at the retina, or higher up in the 
nervous system or brain? 

To address this question, we shall assume case (I) discussed in Sec. 7, that is tq — 
10~^cm and A = 4 x 10~^°s~^. We first ask whether the cis-trans conformational change of 
an individual rhodopsin molecule on absorption of a photon can give reduction with these 
parameter values. A rhodopsin molecule has a molecular weight of about 4 x 10^ nucleons, 
and a diameter of about 4 x 10~^cm. Assuming that in the conformational change all of the 
nucleons in the molecule move by a molecular diameter (likely a considerable overestimate) 
we find from Eq. (9), with n = 4 x 10^, = 1, and £/(2rc) = 2 x 10"^, a reduction 
rate of ~ 3 x 10~^s~^. In other words, reduction with these parameters requires of order 
3 X lO^s, while the conformational change occurs in 200fs = 2 x 10~^^s, so even with the 
enhanced A parameter, reduction cannot occur during the conformational change. For a 
detailed discussion of the possibility of reduction during conformational change, and many 
useful references, see Thaheld [40]. 

On the other hand, with the enhanced A parameter, there is no apparent difficulty 

36 



in reduction occurring during the amplification chain in a rod cell, without invoking signal 
transport in the nervous system to which the rod is linked. According to the review of Rieke 
and Baylor [41], a single catalytically active rhodopsin leads to the closure of several hundred 
ion channels, which over the response time of 300ms blocks the entrance of about 3000 cations 
(sodium or potassium) to the outer segment of the rhodopsin molecule. Assuming that the 
relevant £ here is of the order of the rod cell diameter of 2 microns > rc, we can use the 
estimate of Eq. (8), with N ^ 1 and n = 300 x 3000 x 23 ~ 2 x 10^, giving a reduction 
rate of 2 x 10^s~^. Even if our assumptions here about n and i are optimistic, this estimate 
indicates that with the enhanced A parameter, there is plenty of latitude for reduction to 
be complete by the end of the amplification chain in the individual rod cell. Note that with 
the original CSL A parameter value, the reduction rate would be too small by about three 
orders of magnitude for reduction to occur in a rod cell, requiring the invocation of signal 
transport in the nervous system, as discussed in ref [1]. 

Appendix D: Reduction cirising from phonon emission and thermal fluctuations 

In this appendix we estimate the reduction rate associated with phonon emission by 
electrons and ions, and also the reduction rate associated with thermal lattice fluctuations. 

To start, we note that phonons of wavelength shorter than rc lead to no net center 
of mass displacement of a block of emulsion of linear dimensions rc, and so are not the main 
phonon contributors to reduction. On the other hand, long wavelength phonons can lead to 
a displacement of the center of mass of such a block, and so can potentially have a significant 
effect. Thus, we are interested in making an estimate of the reduction effect of the emission 
of long wavelength acoustic phonons. The effect of short wavelength optical phonons will be 
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taken into account afterwards when we estimate the reduction rate associated with thermal 
lattice fluctuations. 

Consider, then, a lattice block of total mass M with atom site coordinates Ui. For 
a phonon of angular frequency lu, each site will oscillate as Ui = acos{u!t + Si), with 5i a 
site-dependent relative phase, and so the time-averaged velocity is (itf) = c^iiP' jl = u'^i'^, 
independent of i. Since the phonon energy is half kinetic and half potential, we thus have 

-Mcu^e^ = -hcu , (Dla) 

giving for the mean square atomic displacement i'^, 

Muj ^ ' 

For a block composed of N sub- blocks of dimension of the correlation length rc, each con- 
taining n nucleons, we have M — NnrriN, while the reduction rate is given by Eq. (9) in 
terms of N,n,£. Substituting Eq. (Dlb) we get Tr = \f{uj), with 

^ Tin , , 

Note that the factor has dropped out, because when the lattice contains many coherent 
length sized blocks, the N in Eq. (9) cancels against the 1/N imphcit in Eq. (Dlb), which 
reflects the fact that the energy of the single phonon is dispersed over over the whole lattice, 
and so the mean squared amphtude of oscillation of each atomic site scales as 1/N. Note 
also that Eq. (D2a) still assumes that the phonon wavelength is longer than rc, so that all 
nucleons within a block of size rc move together, and the factor in Eq. (9) is applicable. 
To account for the loss of coherence when the phonon wavelength is shorter than rc, a factor 

G{uj) — min[l, (A/rc )^] = min[l, (27rc,/(rca;))'] {D2b) 
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must be included in Eq. (D2a), with Cg the acoustic phonon velocity (the velocity of sound, 
typically ~ 3 x lO^cms"^). This factor takes the value of unity for A > rc, and approaches 
1/n as A approaches the lattice spacing ~ 10~^cm. We will proceed by first doing our 
estimates without including G, and then noting the correction factor, denoted by (G) , when 
G of Eq. (D2b) is included in the integrands. 

Consider now an electron with initial wave number ko that slows down through multi- 
ple emissions of acoustic phonons, as described in Kittel [42] . We wish to integrate Eq. (D2a) 
over the spectrum of emitted phonons, to get the total reduction rate A/tot attributed to the 
emitted phonons. Letting R(^k{t), q) be the rate of production of phonons of wave number q 
when the electron wave number is k{t), we have 

/tot- / dt dqR{k{t),q)f{uj{q)) . {D3a) 

Jo Jo 

The upper limit T of the time integration is determined by k{T) — kth, where kth — 
{3m*kBT)2/h is mean thermal electron wave number (with m* the electron effective mass), 
since once the electron has been thermalized by phonon emission, it is in an equilibrium 
where phonon absorption by the electron is as important as phonon emission. We shall 
assume that both ko and kth are much larger than the threshold momentum for emitting 
acoustic phonons, given in terms of m* and Cg by /cmin = m*Cs/fi. With this assumption 
(which we shall verify) , the upper limit gmax giving the maximum phonon wave number that 
can be emitted by an electron of wave number k is given by gmax = 2 A;. Also with this 
assumption, the time evolution of the electron wave number k is given by dk/dt = —ak^, 
with (7 a constant that will end up dropping out of the calculation. This allows us to change 

39 



variables from t to k, and so Eq. (D3a) becomes 




dka k 



/ dqR{k, q)f{csq) 



{D3b) 



where in the argument of / we have substituted 00 = CsQ. 



The analysis of Kittel [42] also gives an expression for R{k, q), 



R{k,q) 




with a the same constant that appears in the time evolution of the electron wave number. 
So substituting Eq. (D4) into Eq. (D3b), and doing the integrals, we get 



note that this result is essentially independent of the initial electron wave number /cq, as 
long as this is significantly greater than the thermal wave number /cth- Taking the electron 
effective mass m* equal to the electron mass, we get (at room temperature) /cth = lO^cm"^, 
while for the minimum wave number for phonon emission we get /Cmin = 0.3 x lO^cm"^, 
which is much smaller than /cth, as assumed. Evaluating Eq. (D5a) with this value of /cth, 
with rc — 10~^cm and n = 10^, we get /tot — 0.3 x 10^, giving Tr = A/tot ~ 10~^^s~^, which 
is a factor of 10^ smaller than the latent image formation estimate of Eq. (8). This is still 
an overestimate, since we have not included the factor G of Eq. (D2b); when the integrals 
are redone with the factor of G included, we find an additional factor in Eq. (D5a) of 



Thus the reduction rate induced by phonon emission from an electron emitted from a silver 
halide molecule is more than a factor of 10^ smaller than Eq. (8). 
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{Dhh) 



The silver and bromine ions that diffuse to the surface of the grain will already 
be in thermal equilibrium. We assume that they obey a dispersion relation of the form 
E — p^/(2m*), with m* an effective mass similar to the ionic mass, and so their mean wave 
numbers will be given in terms of this effective mass by kth. = {3m*kBT)^ /h. For a particle 
of effective mass m* the ratio /Cmin/^th scales as (m*)^, and so for an ion of atomic weight 
~ 100 we get by comparison with the electron case calculated above, 

A;min/A;th ~ (2 x 10^)5 x 3 x 10"^ ~ 10 . {D6) 

Hence only the extreme tail of the thermal velocity distribution for these ions is above 
threshold for phonon production, leading to a suppression of the associated reduction rate 
by a factor exp{-1.5k'^^^/ k^^) ~ 10~^^. 

When both electrons and interstitial ions are in thermal equilibrium, they continue 
to emit and absorb phonons. The associated reduction rate will be driven by the thermal 
average of ^^iuj — lif)^, with u]'^ the atomic coordinates in superimposed states in the 
wave function labeled by the respective indices 1,2. An upper bound on this sum should 
be given by ignoring expected cancellations in the individual terms (that is, for most atoms 
we expect u] ~ u'^), and replacing the sum by twice the thermal average of X^j^i- The 
associated reduction rate is given in turn by the twice the integral of Eq. (D2a), including 
the factor G of Eq. (D2b), over the Debye phonon spectrum. When G — 1, the result is 
Fr = 2 A/tot, where 

_9 ksTn^N 

with ud the Debye frequency. [A similar result, with ujj replaced by Scu^p^, with Uopt the 
lattice vibration frequency, is readily obtained from Eq. (9), when £^ is taken as the thermal 
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average of the atomic position, given by (u^) — 3/cBT'/(Matom<^opt)-] Including the coherence 
factor G in the integral over the Debye spectrum gives an extra factor 

(G) ^Sncs/ircujo) ■ {D7b) 

Evaluating Eqs. (D7a) and (D7b) at room temperature, with rc — lO^^cm, n — 10^, Matom = 
lOOmiv, and a;D = 3 x lO^^s"^, we get 

Tj? ~ 2 X lO-^TVs"^ , {D8a) 

or with N = 20, 

Tfl - 4 X IQ-^s"^ , {D8b) 

a factor of three bigger than the latent image formation estimate of Eq. (8). As noted above, 
Eq. (D8b) is expected to give a substantial overestimate of the contribution of thermal 
fluctuations to the reduction rate, since for most atoms one expects (u] — uff « {ulY + 
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